The problem of phase breaking in the electronic conduction in mesoscopic systems: a 

linear-response theory approach 
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We study the problem of electronic conduction in mesoscopic systems when the electrons are 
allowed to interact not only with static impurities, but also with a scatterer (a phase breaker(PB)) 
that possesses internal degrees of freedom. We first analyze the role of the PB in reducing the 
coherent interference effects in a one-electron quantum mechanical system. In the many-electron 
system we can make a number of quite general statements within the framework of linear-response 
theory and the random-phase approximation. We cannot calculate the conductivity tensor in full 
generality: we thus resort to a model, in which that tensor can be expressed entirely in a single- 
electron picture. The resulting zero-temperature conductance can be written in terms of the total 
transmission coefficient at the Fermi energy, containing an additional trace over the states of the 
PB. 



I. INTRODUCTION 

The scattering approach to quantum electronic trans=. 
port in mesoscopic systems was devised by LandauerEJ 
and later extended by a number of other authors (see, 
for instance, as representative articles, Refs. ||J| and the 
references contained therein) . In an independent-electron 
picture, it aims at understanding the electric conductance 
of a sample in terms of its scattering properties. The 
problem of electric transport is thus converted into that 
of solving the quantum-mechanical scattering problem of 
an electron that impinges on the sample through leads 
that, ideally, extend to infinity, once the experimental 
environment the sample is attached to in the laboratory 
is disconnected. An approach to this problem using the 
methods of linear- response theory (LRT) has been given, 
for instance, by the authors of Refs. ( see a l so other 
publications referred to there). 

In the original conception of the scattering approach 
to electronic transport, inelastic electron scattering or 
other phase-breaking mechanisms are not allowed inside 
the sample. As a result, the phase of the wave function 
is completely coherent in that region. Yet, in various 
circumstances the effects of the electron-electron interac- 
tion or the interaction with the phonon field may apt be 
negligible. In a further development of the theorya, the 
single-electron picture is maintained and phase-breaking 
events in a given region are modelled by attaching the 
system to a "fake wire" , which in turn is connected to 
a phase-randomizing reservoir, so that there is no phase 
coherence in the wave function for an electron entering 
and exiting the reservoir. The chemical potential of the 
reservoir is chosen so that the net current along the fake 
wire vanishes. This model provides sensible answers and 
has been used, for instance, in, the study of electric trans- 
port through quantum dotsOil 

A number of authors have attempted to generalize the 
scattering approach to include inelastic scattering explic- 



itly, instead of modelling it as described above. In Ref. 
|8| the problem of quantum transport in the presence of 
phase-breaking scattering is formulated using an exacly 
soluble model for the electron-phonon interaction and, 
using linear-response theory, a generalized conduction 
formula is found. Ref. ^ uses Landauer's approach and 
analyzes the effect of a single impurity scatterer, at which 
both elastic and inelastic processes can occur. The form 
of the electron-impurity interaction can be quite general; 
in the weak-scattering limit, in which Born approxima- 
tion is invoked, the authors arrive at a generalized con- 
duction formula. The possibility of energy exchange with 
the scatterer makes Pauli blocking effects important and 
an extension beyond Born approximation seems difficult. 

The authors of Ref. [ll^ analyze in great generality the 
problem of quantum-mechanical phase breaking: they 
consider the interference between the terms of an electron 
wave function arising from two different electron paths 
and study the effect on that interference of an "environ- 
ment" the electron can interact with. In an experiment 
in which the electron and not the environment is mea- 
sured, the coordinate of the latter is integrated upon: as 
a result, the interference is lost if, in the two interfering 
partial waves, the states of the environment are orthogo- 
nal to each other. It is emphasized that, for this to occur, 
energy exchange between the electron and the environ- 
ment need not be invoked. That this loss of interference 
is irretrievable is a quantum- mechanical effect, common 
to a number of situations, like the two-slit experiment. 

In the present paper we plan to incorporate the mech- 
anism of Ref. [l0| -briefly described in the previous 
paragraph- to the LRT approach to transport provided 
by Refs. |||| in order to study the problem of phase 
breaking in the electronic conduction in mesoscopic sys- 
tems. The emphasis on the mechanism of Ref. [lO] for 
phase breaking is the main difference between the present 
and previous work on the problem. Another characteris- 
tic of the present paper is that the applications are dis- 
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cussed in terms of the scattering matrices of the various 
impurities (static or non-static), a procedure that has 
been found adva ntageo us in a number of previous publi- 
cations (see Refs. [11 12 and references contained therein). 

Quite generally, we can pose the problem by think- 
ing that the electrons, besides suffering elastic collisions 
with static scatterers, interact with a number of scatter- 
ers, or phase breakers (PB), that possess internal degrees 
of freedom and can live, say, in m possible quantum- 
mechanical states altogether. Even in the absence of 
the electron-electron Coulomb interaction, the problem 
is now no longer one of a single-electron, but is a full 
many-body problem: one electron is incident on the PB, 
this being in some state /j,; after the interaction there is a 
certain probability to find the PB in state v, and this is 
what the next electron coming in will see. This memory 
effect, or, equivalently, the electron-electron interaction 
induced by the PB, gives rise to_ a, situation similar to 
the one found in Kondo problenxj and we are bound to 
find similar complications. 

The paper is organized as follows. In order to become 
acquainted with the physical phenomena produced by a 
PB we first study, in the next section, the problem of a 
single electron scattered by two static impurities and a 
PB. We show how the interference terms that occur with 
the static impurities alone are affected by the presence of 
the PB. The discussion parallels that given in Ref. [To|and 
is done in terms of the scattering matrices of the various 
impurities. 

In Sec. Ill we pose the conduction problem (a many- 
electron problem) of an electronic system with static im- 
purities and a PB [but not subject to a magnetic field, 
so that we have time-reversal invariance (TRI)], from the 
standpoint of LRT. We show that we can make a number 
of quite general statements. However, we reach a point 
where we are unable to calculate the conductivity tensor 
in full generality: we thus resort to a simplified, soluble 
model that we introduce in Sec 



Our conclusions are discussed in Sec. VI 



II. QUANTUM INTERFERENCE IN A 
ONE-ELECTRON SCATTERING PROBLEM: 
THE EFFECT OF A PHASE BREAKER. 

We analyze in this section the scattering problem of 
a single electron interacting with a combination of three 
scatterers in series: two static ones and a PB in the mid- 
dle. We show how the interference terms that would 
occur with the static scatterers alone are affected by the 
presence of the PB. 

For simplicity, the problem is treated as a ID one. It 
is described by the Hamiltonian 



2m 



- Vi{x) + J2 lM> V^{x) H + V 3 (x). (2.1) 
1.1. i j 



In this equation, V\ (x) and V3 (x) are the potentials aris- 
ing from the two static scatterers, and the third term 
represents the interaction of the electron with the PB; 
the m states of the latter, denoted by /i, v, are degener- 
ate in energy. Below, we shall find it convenient to write 
our states as m-component "spinors" , so that H acquires 
a matrix form, where rows and columns are associated 
with the m states of the PB. 

Instead of modelling the potentials for the various scat- 
terers, we have found it advantageous to model their scat- 
tering properties through the corresponding scattering 
matrices. In this language, the three scatterers are de- 
scribed by the S matrices Si, S2 and S3, respectiverly. 
We write the S matrix S2 for the PB as 



S 2 = 



T2 t' 



(2.2) 



tensor is calculated within that model and is found to be 
expressible entirely in terms of a single- electron picture, 
i.e. in terms of single-electron Green's functions. It is 
then shown that the resulting zero-temperature dc con- 
ductance can be expressed in terms of a total transmission 
coefficient at the Fermi energy ep, but no w con taining a 
trace over the m states of the PB, Eq. (4.34). Within 
the restrictions of the model, the result does not depend 
on the strength of the electron-PB interaction, as other 
analysis do. We present (see discussion at the end of 



Fv. The conductivity where the reflection and transmission matrices (for inci- 
dence from the left or from the right, respectively) r2, r' 2 , 
t2, t' 2 are m-dimensional, with matrix elements r% , etc. 
The S matrices for the elastic scatterers on the left and 
on the right of the PB are written, respectively, as 



Si 



s. 



section III around equation (3.17)) a speculation as to 
the validity of our main result beyond the assumptions 
of the soluble model, in the strict linear-transport regime 
and above the Kondo temperature (which is taken to be 
extremely low) associated with the m-level PB. 

In Sec. we set up a random-matrix description of the 
electron-PB system, with possible applications to chaotic 
cavities, in order to calculate the effect of the PB on the 
average conductance and its fluctuations. The limita- 
tions of the model become apparent here. 



rih 
til„ 



r 3 I„ 



¥ T 
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l '3 J m 



1 3 ^"3 



(2.3) 



(2.4) 



where n, ... , r%, ... , are just complex numbers (we are 
in ID) and I m is the m-dimensional unit matrix in the 
space of the PB states: recall that scatterers 1 and 3 do 
not change the state of the PB. 

The total transmission matrix t for the chain of three 
scatterers in series is given by 



t = (t 3 I m ) ■ 



-h 



r 12 r 3 



(2.5) 
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In this equation, ti% is the transmission matrix for the 
system formed by S\ and 5*2, given by 



ti2 — t 2 



1 



(r£)r 2 



{til m ) 



(2.6) 



Similarly, r' 12 , the reflection matrix for the combination 
Si, S 2 , is given by 

1 



1 m V l) '2 

To be more specific, we now choose r 2 



(2.7) 
r' 2 = in 

the matrix S 2 of Eq. (|2j) that defines the PB. Thus the 
PB will not change the electron momentum; however, we 
shall see that it reduces, in the electronic current, the 
interference among the multiply reflected paths occurring 
between the two elastic scatterers. As a consequence of 
this choice, the matrices t 2 and t' 2 are m-dimensional uni- 
tary matrices 



t 2 t 2 — I m , 



Also, ti2 reduces to 



and r' 12 to 



t 2 t 2 



tati 



r'12 = t2r[t' 2 



(2.8) 
(2.9) 
(2.10) 



We thus have, for t 

t = (t 3 ti) ■ 



{hh)t 2 



i 



(r[r 3 )t 2 t 



(*3*l)*2 



T t 2 



1 



' I m - (r[r 3 ) t' 2 t 2 "* I m - an 

where we have defined the complex number 
a = r[r 3 = \a\ exp(ip) 



and the unitary matrix 



u = t' 2 t 2 . 



(2.12) 



(2.13) 



The total transmission matrix t of Eq. (2.11) is Tri- 
dimensional. Its element t^ v gives the probability am- 
plitude for the process: {the electron comes from the 
left, the PB being in state v\ — > {the electron is trans- 
mitted to the right, the PB being shifted to the state 
/i}. The corresponding probability is T^ v = \t^\ 2 . Now, 
T" = ^2 T tiv is the transmission probability when the 
PB is initially in state v and the PB is not observed; T v 
can be written as 

t v = {thy v 



where Eq. (2.8) was used. Let us emphasize that the 
superscript v specifies the initial pure state of the PB. If, 
however, the PB is initially in a mixed state (a situation 
of pa rticular interest for the conductance problem, Sec. 
IV 1), then an additional sum over v is needed. Assum- 
ing that the PB can be found, with equal probability, 
in each of its m states, one obtains for the transmission 
coefficient T: 



r = -Vr = —trtth) 



= TiT 3 —tr 



(2.15) 



(2.16) 



I m - a*ut I m - a 
It is instructive to expand t of Eq. ( |2.11 ) as the power 



t=(t 3 ti)t 2 [I„ 



2 2 

+ au + a u 



(2.17) 



This series can be easily interpreted in terms of the multi- 
ple scattering processes that occur between the two elas- 
tic scatterers, influenced by the PB in the middle. The 
transmitted wave function is a linear combination of all 
these terms, arising from internal multiple reflections. 
We now study a number of particular cases. 



We set 



A. The Case m=l 



t 2 = t'= 1. 



, (2.11) From Eq. (2.11) we have 



*3*1 



1 



r[r 3 



Ui 
l-a 



(2.18) 



(2.19) 



In this case, scatterer 2 is not a PB, but a static scat- 
terer: it is thus like having just the two elastic scatterers 
[more generally, we could choose t 2 = e m ,t' 2 = e 1 ' 3 ; then 
the effect of scatterer 2 would simply be the addition 
of the relative phase (a + (3) between the original scat- 
terers 1 and 3; that extra phase could equally well be 
obtained, for instance, by setting the two elastic scatter- 
ers a distance d farther apart, where kd = (a + /?)]. The 
transmission probability discussed above is 



T 1 =T 



T 3 Ti 



= T 3 Ti 



1 



l-a 2 



^ m 1 + \a\ A + 2Rea _ m 

= h aia = lcoh - 

1 - a 2 



(2.20) 



T3T1 



I m - a*ut I - , 



(2.14) 



This result will be termed the fully coherent response 

T C oh ■ 
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B. The Case m=2 

Now the PB has two (orthogonal) states. We follow 
the various multiply scattered terms occurrin g bet ween 
the two elastic scatterers -as given by Eq. (2.17)- in 
order to understand more closely what the PB does to 
them. We consider two examples. 

1. Let 



to 



t'o 



so that, from Eq. ( |2.13| ) 

u = t' 2 t 2 = h 



1 

1 



1 
1 



(2.21) 



(2.22) 



We examine the multiple-scattering series (2.17). In each 



passage through the PB, the state of the latter is shifted 
to the orthogonal state. But, after the pair of reflections 
described by the product a = r[r 3 , the PB is visited twice 
and is then back to the original state. In other words, in 
each passage, the PB exactly undoes what it did in the 
previ ous o ne. This is the significance of u = t' 2 t2 = li in 
Eq. (|2~22|). 



We thus find for t, Eq. 



t = (t 3 ti)t 2 



1 / ' 

1 - 7-^3 



(2.23) 



which leads to Eq. ( 2.20| ) , exactly as for the case with no 
PB. 
2. Let 



2 l 2 



1 

1 



(2.24) 



This could be obtained, for instance, with the choices 



or 



*2 = to = - 



1 + i l-i 
1 - i 1 + i 



(2.25) 



(2.26) 



Now u shifts the two PB states: i.e., state |1) is shifted to 
|2) and viceversa. This fact has imp ortant consequences. 
The multiple-scattering series (2.17) for t now gives 



t = (t 3 h)t 2 



a 2 h 



3 

a a a 



(2.27) 



which divides naturally into an even-order and an odd- 
order contribution, which can be summed up to give 



t=(t 3 ti)t 2 



1 



I -a 1 



1-a 2 



(2.28) 



Suppose that in the incident state the electron comes 
from the left and the PB is in the pure state 



|0> 



(2.29) 



say. The transmitted wave function on the right is thus 



*3*1 _ e ikx 



2(1 -a 2 ) 



(l + a)(a + /3) + (l-a)(a-/3)i 
(1 + a)(a + /3) - (1 -a)(a- 0)i 



(2.30) 



The transmission probability, that we call T , is given by 



l + |a| 2 + 2Re(a/3*)2Re(a) 
II -a 2 ! 2 ' 



(2.31) 



In the absence of the PB, on the othe r ha nd, we have 
the fully coherent response T C oft_of Eq. ( |2.20 ). The inter- 
ference term 2Re(a) in Eq. (2.20) has been reduced, in 
Eq. (2.31 ), by a factor, whose magnitude |2Re(a/3*)| < 1. 



The effect of the PB, and of having measured the electron 
but not the PB, has thus been to decrease the magnitude 
of the interference term. We stress again that this ef- 
fect i s there even if the incident state is the pure state 
( p. 29 ). In particular, for either pure state a = 1, ft = 0, 
or a = 0, P = 1, we obtain, in the notation introduced 



right after Eq. (2.13) 



Jll rj-iZ 



r 3 Ti- 



II -a 2 



(2.32) 



where the interference term in the numerator is absent. 
This last result could be obtained directly from Eq. 
( |2.14 ) . For a mixture of these two states we thus obtain 
the same answer, i.e. T = T 1 = T 2 . 



C. The Case m — > 00 

Let us inquire as to what kind of a PB would lead to the 
"classical" composition rule for the two elastic scatterers, 
i.e. to the equation 



T = 



T1T3 



1 — R\R 3 



(2.33) 



We shall see that for this to occur we need a PB with 
many states, i.e. m — > 00. We analyze two possibilities. 
1. Choose 



(t2 



(*2 



2 //if 



(2.34) 



For finite m we use periodic boundary conditions, i.e. 
r?i+l —> 1. But below we are interested in m — » 00. 
Assume that the electron is impinging from the left and 
the PB is initially, say, in its first state |1). Each passage 
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of the electron through the PB will flip the latter to its 
next (orthogonal) state. Therefore, from the expansion 
( [2.17 ), the transmitted wavefunction is: 



l*fc)l 



e ifca: t 3 ii[|2) + a|4> + a 2 |6) 



(2.35) 



It is now clear that the transmission coefficient T 1 for 
t he el ectron (without observing the PB) is given by Eq. 

(HI). 



2. We now discuss a model that allows varying the 
degree of dephasing in a continuous fashion, thus per- 
mitting the description of the crossover between the fully 
coherent r espon se ( [2.20 ) and the fully incoherent, or clas- 
sical, one ( 2.33 ). 

We choose 1% — t' 2 — exp(iH/2) where H is an m x m 
Hermitean matrix. Its eigenvalues 9 are chosen so that 
their density mg(9) has width A (a trivial way to get that 
is to take a diagonal matrix and put the required spec- 
trum by hand, then make an arbitrary unitary transfor- 
mation). Obviously, the eigenvalues of t\ will be exp(i9), 
where 9 is an eigenvalue of H . Assumi ng th at the PB is 
in a mixed state and, thus, using Eq. (2.15), we have 



T = T X T Z 



1 



\ " 



1 



TO ^— ' II 

r=l 1 



ae' 



(2.36) 



which in the m — > oo limit can be written as 

g(9)d9 



T = TiTg 



o |l-ae*| a 



(2.37) 



Rather then attempt to do the last integral exactly for 
a, give n g(9), it is very instructive to expand, as in Eq. 
( 2.17 ), each of the two factors in the denominator as a 
sum of powers such as a k exp(ik9), and the complex con- 

k' 

jugate, (a*) exp(— ik'6). All th e "diagonal" terms give 



the classical result, Eq. ( p. 33 ), as above. The mixed 
terms with k ^ k! give the interference. With vanishing 
A, we get the full interference terms. Otherwise the in- 
terference terms are killed once \k — k'\X > 0(1). Thus 
1/A plays the role of the dephasing length, in units of the 
distance between the static scatterers. 

It is possible, using the above, to analyze a number 
of interesting physical situations. Consider a real and 
close to unity, i.e. a = 1 — 6, with < 6 <C 1 and g(9) 
picked at 9 — 0. Take, for example, the peak of the reso- 
nant transmission for the fully coherent case, A = 0. The 
height of the resonance is ~ 1/6 2 . The relevant number 
of bounces (the number of significant terms in the series 
obtained upon expansion of 1/(1 — ae 10 ) needed to form 
the resonance) is ~ 1/6, Now introduce a finite A. If 
A <C 5, it will affect distant (i.e. further than 1/6) terms 
in the series and thus will have no effect on the height of 
the resonance. Dephasing starts to "hurt" the peak once 
A»<5. 



III. THE CONDUCTION PROBLEM IN THE 
PRESENCE OF A PHASE BREAKER IN 
LINEAR-RESPONSE. 



In this section we discuss the conduction problem -a 
many-electron problem- in the presence of a-EB from the 
standpoint of linear-response theory (LRT)EJ, following 
the random phase approximation (RPA) scheme devel- 
oped in Refs. ^ and pi The system to be studied has the 
geometry shown in Fig. [l|: the constriction represents 
the sample, where we allow for the presence of a PB. As 
usual, the expanding horns represent the external leads 
that, in a laboratory setup, are attached to macroscopic 
bodies. 



dS'- 




- dS' 



FIG. 1. The geometry for the electronic conduction problem 
studied in the text. The constriction represents the sample and 
the expanding horns, the external leads. Beyond the surfaces S+ 
and S— the potential 5^> w takes on constant values. 

When no external voltage is applied, the whole system 
is in equilibrium and is described by the unperturbed 
Hamiltonian 



i 



*<3 



ti-r s \ 



2m ^ Irj — R/l ^ |r; - R, 



(3.1) 



Here, — e is the electronic charge, the position variable 
of the i-th conduction electron, R/ the (static) position 
of the J-th ion with positive charge ej (screened by the 
bound electrons), and ej and Rj the charge and position 
of the j-th impurity, considered to be static. Thus H° 
contains the kinetic energy, the interaction of the elec- 
trons with the ions and the static impurities, as well as 
the electron-electron interaction, for which no approxi- 
mation is assumed for the ti me being. It also contains, 
in the last two terms of Eq. (3.1), the intrinsic Hamilto- 
nian of the PB and its interaction with the electrons, of 
the type introduced in the previous section for one elec- 
tron. We do not consider a static magnetic field to be 
present, so that the problem is time-reversal invariant. 

We denote by po(r) and 4>o(r) the equilibrium charge 
density and potential, respectively, that satisfy Poisson's 
equation 
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V 2 o (r) = -47r j o (r). 



(3.2) 



Application of an external voltage with frequency u> will 
cause a current density J£(r) in the system. It will also 
lead to a change in the charge density and in the po- 
tential. We denote these changes by 8(f(r) and 84^(r) 
and empasize that these are full changes, with respect 
to the equilibrium values po(r) and </>o(r). In the RPA 
approximation, which is employed in the present paper, 
there is no need to separate the full changes into external 
and induced parts. The essence of the RPA is to omit 
the electron-electron interactions from the Hamiltonian 
Hq but to use, instead, the full potential 5<j> u (r) (rather 
than the external one) in the formulation of LRT. The 
full potential is then determined self-consistently from 
the Poissan equation. We thus have the three following 
equations^: 



J-(r)= / d 3 rT-(r,r')^(r'), 



5p u (r) = J d 3 r'n w (r,r')(5^(r'), 

V 2 6^{r) = -4rrV(r) 
with the kernels 



(3.3) 



(3.4) 



(3.5) 



r£(r,r') = -- / dre MT - 7T 
n Jo 



n"(r, r') 



ia(r,r),p(r',0) 



p e l(r,T),p e i(r',0) 



00 

(3.6) 



00 

(3.7) 



where the limit 7 — > is implied. The expectation val- 
ues of the interaction picture operators are taken with 
respect to the Hamiltonian 



Hon — H — 



i<j 



with the electron-electron interaction switched off. 



(3.8) 



The Poisson equation (3.5) should be supplemented by 
boundary conditions. Sufficiently far inside the horns the 
local current density J% (r) becomes vanishangly small, so 
that these distant regions remain practically in equilib- 
rium. This means that well inside the horns 5<f) u (r) ap- 
proaches constant values, o~4> u {— 00) and d^i+oo). The 
difference 



; (-oo) -S(/> u (+oo) = V u 



(3.9) 



is the total potential drop on the system sample+horns 
(V™ generally does not coincide with the external EMF, 



since some voltage drop can occur in other parts of the 
circuit, e.g., near the points where the external EMF 
source is connected to the horns) . At the internal bound- 
ary of the system one should require zero current density 
normal to the boundary 



J n (r s ) = 0, 



(3.10) 



r s being a poin t on th e int ernal boundary. 

Solv ing ( 3.3), (3.4), ( |3.5|) with the boundary conditions 
( |3.9[ ), (3.1C) would enable one to determine the charge, 
current density and field distribution within the system. 
One could then compute the total current I u and, thus, 
the conductance = I u /V u . This is a formidable prob- 
lem. A great simplification, however, occurs in the dc 
limit, thanks to a result obtained in Ref. [if] that states 
that, in that limit and for fixed r, r', the conductivity 
tensor is divergenceless; i.e. 



tyoZp (r,rO-0, 



(3.11) 



The conductivity tensor relates current density to the 
elect ric field [rather than to the potential, as in equation 
©]; i-e. 



J» = - d 3 r'a» (r,r')d> 5cl> u (r'), 



(3.12) 



and is given in terms of the current-current correlation 
function as 



1 



e 2 n (r) 



vmoj 



dre 1 



ja(r,r),^(r',0) 



00 
(3.13) 



no( r ) bei ng th e electron density in equilibrium. Integrat- 
ing Eq. ( 3.12 ) by parts shows that, in the dc limit, the 
current density is insensitive to the full potential profile 
within the sample and depends only on the total poten- 
tial drop between the two distant surfaces well inside the 
horns (Fig. [l|): 



j: 



LJ^O 



(r) = - [<^°(+oo)r+(r) + <5<^°(-oo)r-(r)] , 

(3.14) 



where 



r£(r)= / 



(rX). 



(3.15) 



This observation paves the way for a derivation of a Lan- 
dauer formula from the LRT, and demonstrates that in- 
teractions, within RPA, do not affect the conductance. 
This conclusion, known for purely elastic scatterersU, re- 
mains valid also in presence of the PB. 

For zero temperature, and in the absence of PB, the 
answer is the well known one 



G 



G 



E 



(3.16) 



ah 



where t is the single-particle transmission matrix at the 
Fermi energy, from well inside the left horn to deep inside 
the right one. 

When the system is descr ibed by the unperturbed 
Hamiltonian Hqo of Eq. (3.8), the electrons can change 
the state of the PB through the e-PB interaction rep- 
resented by the last term in Eq. ([O]); the coupling of 
the electrons with the PB has the consequence that the 
former are no longer independent. This results in a com- 
plicated structure for the TV-electron eigenstates of H ao , 
and the calculation of a\ a (r, r') is, in principle, no longer 
feasible along the lines followed in the absence of the PB. 
In the next section we discuss a model for the e-PB in- 
teraction that does lead to a solution along similar lines. 
One might speculate on physical grounds that our final 
result, equation (4.34), should be valid under the follow- 
ing assumptions. One should first take the temperature 
to be very low, but much larger than the Kondo tem- 
perature, Tk, due to the interaction of the PB with the 
electron gas. Furthermore, one should stay in the strict 
linear response regime, where for a finite distance L be- 
tween the two reservoirs, the current satisfies: 



e/I > Ljvp. 



(3.17) 



This assures that the separation in time between consec- 
utive electrons participating in the transport is so large 
that the first electron reaches the downstream reservoir 
and thermalizes there before the next electron starts its 
journey. This can be expected to eliminate the electron- 
electron interaction mediated by the PB, which is a co- 
herent second-order process. 



IV. A SOLUBLE MODEL 



or, in terms of its matrix elements 

^2 



*8Z = E 



(4.3) 



We choose A' 7 constant, i.e. independent of a (and hence 
we set it equal to zero) , so as not to have constant terms 
in the off-diagonal matrix elements (i / c. In the lan- 
guage of Eq. (3.1), the energies of the PB states are 
degenerate and set equal to zero. 

The Schrodinger equation in the D and ND-basis is 



ifoo |*) = E\V), 



(4.4) 
(4.5) 



respectively, with 

|*) = 0|*>. (4.6) 
In the D-basis, the Hamiltonian Hqq can be written as 



with 



and 



H oo W = Ho (i)6 a 



In the ND basis, 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



We assume that Hqo of Eq. (3.8) has such a structure 
that, in a suitable PB basis, labelled by a, a' below, it 
acquires the diagonal form 



H 



00 



E 



(4.1) 



Just as we did in Sec. f| we write the Hamiltonian in 
matrix form, where rows and columns, i.e. cr, a', label 

the to states of the PB. The bar in H 00 indicates that 
the Hamiltonian is expressed in the new PB basis, that 
we shall call, for short, the D-basis (D for diagonal), as 
opposed to the original, or ND-basis. 

In the ND basis, H 00 is obtained from H o by means of 
a constant, real (in order to preserve reality of the Hamil- 
tonian, and hence time-reversal invariance) orthogonal 
transformation O (an m-dimensional matrix), i.e. 



Hnn — OHr,nO T , 



(4.2) 



with 



and 



H, 



u^(n) = Vo" ff c/ ff (r,)o'- 



cr=l 



(4.11) 



(4.12) 



If the matrix = [O^f has nonzero determinant, 
we can find the t/ <7 (ri)'s (a =1,- • •, m) that reproduce 
any given set of diagonal potentials J7 MM (ri)'s (fj, =1,- • ■, 
to); but then we are left with no freedom to select the 
off-diagonal interactions {/'"'(rj) (/i ^ v), which become 
uniquely determined by the i/' T (rj)'s. So, i t is clear that 
the matrix elements C/ MI/ (ri) of Eq. (4.12) show strong 
correlations among themselves. These correlations make 
it possible to find a D-ba si s in which the Hamiltonian 
takes the form of Eqs. (4.7)-(4.9) and the problem breaks 
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up into m independent single-particle ones: this is the 
feature that makes the problem soluble. Solving the 
Schrodinger equation (4.4) thus reduces to solving the 
m single-particle Schrodinger equations 



^(r,) = «(r i ), 



a=l,- 



(4.13) 



We shall assume that none of the £/ CT (ri)'s admits bound 
states. 

In the D-basis, the states 
















a = 1, 



(4.14) 



with a nonzero value in the er-th entry, form a com- 
plete set of orthonormalized (in a (5-function sense) 
single-particle states, eige nfunctions of the single-particle 
Hamiltonian matrix (4.S). 

In the D-basis, the S matrix associated with a scatter- 
ing solution has the form 



s =s s„ 



(4.15) 



If the problem admits N open spatial channels, each S a 
is 27V-dimcnsional. In the ND-basis S takes the form 



(4.16) 



and is 27n7V-dimensional. The number of independent 
parameters associated with each unitary symmetric ma- 
trix ST is N(2N + 1) and is thus mN(2N + 1) for the 
total S and hence for S. This makes it clear that the S 
matrices allowed by our soluble model do not have the 
"generic" structure used in some of the considerations 
of Sec. ||, but have a rather restricted one: in fact, a 
generic miV-dimensional S-matrix has a larger number, 
i.e. mN(2mN + 1), of independent parameters. 

For two particles, say, we have the states (not antisym- 
metrized yet), with a = ±1 







ki k.2 



o 





(4.17) 



To antisymmetrize we use a second-quantization lan- 
guage, so that for N electrons we have the states 









o 







|o>, 



(4.18) 



with (c£)t creating one electron in state ^(r) and |0) 
being the electron vacuum. 



1. The conductivity tensor and the conductance 

The conductivity tensor cr^ (r, r') is given in Eq. 
( |3.13 ). The expectation value occurring in that equation 
has to be understood as 



ja(r,T),jp(r',Q) 



00 



= J2 P(MM;a)(AfM;a [j a (r,T),j (r',O) 

NMcr 



ATM; a 
(4.19) 



Here, \Af; a) are the states of Eq. ( [4.18| ) in the D PB 
basis, Af being the number of electrons and M an abbre- 
viation for the configuration fej., &2, • ■ •, ^V- The states 
\ AfM ; <j) are eigenstates of the Hamiltonian Hqq of Eq. 
Q, with the energy E NM ; a =4 1 +--- + ^k N - 

The current operator j a (r,r) in the interaction repre- 
sentation is given by 



j a (r,r) =e* HoaT j a {r)e-* 



-H 00 T 



(4.20) 



In the D PB basis t he H a mil tonian ffoo has the diagonal 
form given by Eq. ( |4.7| )-( f4~9| ); since the current operator 
j a {r) does not depend on the PB explicitly, j a (r, t) takes 
the diagonal form 



where 



fa (V,T) = %(T,T)S a 



jZ{t,t) = e« H °° T j a (r)e s- H °° T , 



(4.21) 



(4.22) 



H 00 being gi ven b y Eq. ( |4.7| ). 

Definition (4.19) implies, as usual, that the density ma- 
trix is diagonal in the representation in which the Hamil- 
tonian is diagonal, with diagonal elements P(AfM;a). 
Explicitly, P(AfM; a) is given by the grand-canonical en- 
semble (understanding now the labels TV, M as the set 
of occupation numbers n\, n^, ■ ■ ■ ) as 
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P{n 1 ,n 2 , ...;cr) 



-/JEHU 



Z(/3, M ) 

where the grand partition function is 

m 



ith 



(4.23) 



(4.24) 



(4.25) 



Wc define the conditional occupation probability 



P(n 1 ,n 2 , ...|cr) = 



(4.26) 



the condition being that the PB be precisely in the state 
a of the D basis; its trace is 1. Then 



P(m,n 2 , a) = p(a)P(n 1 ,n 2 , ...\a), 



where 



p(a) 



Z(P,n;v) 



(4.27) 



(4.28) 



is the probability of finding the PB in state a, the tem- 
perature and chemical potential, not indicated explicitly, 
being /3, fi. Of course we have 



(4.29) 



and so the trace of ( 4.27 ) is 1 



We can thus write the expectation value (4.19) as 
[7«(r,r)J /J (r',0)]) oo 

rn 

= ^2p(v) P(n u n 2 ,...\a) 

(7—1 Hi,Tl.2)"' 

; ni,n 2 ,...;cr|^(r,r),^(r',0) 



nx,n 2 , ...;<r 

(4.30) 



Thus in the D PB basis the problem breaks up into 
m independent problems, for the Hamiltonians H 00 , 
er = 1, • • •, m. The conductivity tensor a^a (r, r') of Eq. 
( |3 .13 ) can thus be written as 



^ (r,r') = ^p(aK^(r,r'). 



(4.31) 



where (r, r') is a conductivity tensor that can be 
expressed-,|«jitirely in terms of single-electron Green's 
functionsEnia for the Hamiltonian H 00 (i). 



We write the conductance G in terms of the "dimen- 
sionless conductance" g as 



(4.32) 



As to — ► and then the temperature — ► 0, we find, for 
"spinless electrons" 



g = T = TjTpi*)tr[ff')ir], 



(4.33) 



(7 = 1 



t a being the transmission matrix (an N x N block of 



the matrix S of Eq. ( 4.15 )) arising from the potential 
[/ CT (r); the trace in the above equation is over spatial 
channels, as usual. Should U a (r) = U(r), i.e. indepen- 
dent of a (and hence U >iv {y) = U{y)8 i11 , in any basis), 
the above formula would go over into the standard one. 
In the ND PB basis we finally find (a, b being spatial- 
channel indices) 



fivis' ab 



- V- 
PPB 



Tr (ppBtU) , 



(4.34) 



the trace being now over the spatial channels and the PB 
states. Wc have defined 



(4.35) 



Eq. ( 4.34 ) has the structure of the standard Landaucr 
formula, with an extra avera ge ov er the PB states. This 
is also the structure of Eq. (2.15), that was obtained in 
the study of a single electron interacting with a PB with 
equal weights assigned to every PB state. 

Should there be circumstances where the various p(cr) 
discussed above, at zero temperature, were all equal to 
l/m, we could write 



/i,i/=l ab 



(4.36) 



We recall that i t is f or an e-PB interaction of the type 
described in Eq. (4.12) that 



1) the many- electron Hamiltonian (with the e-e inter- 
action switched off) breaks up, in the D basis, into m 
independent single- electron Hamiltonians, and hence 

2) the conduction problem breaks up in a similar man- 
ner and is soluble in terms of single- electron quantities. 

We also remind the reader that within our model the 
single-e-PB S matrix, and henc e the transmission am- 
plitudes appearing in Eq. (4.34), are not "generic", but 
have a rather restricted structure. 



The final answer, though, i.e. Eq. (4.34), is a very 
appealing one and is likely to be valid beyond the situa- 
tion envisaged by the present model, i.e. to cases where 
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the e-PB S matrix has a more general structure. How- 
ever, in such a generic case it is not known to us un- 
der what approximations the conduction problem -even 
with the e-e Coulomb interaction switched off- could still 
be reduced to independent single-electron problems and 
thus expressed in terms of single-electron quantities. The 
best we ca n do at present is conjecture the validity of 
Eq. ( 4.34 ) for a generic e-PB S matrix, under suitable 
approximations, as ex plained at the end of section III 
around equation (3.17). 



V. A RANDOM-MATRIX MODEL FOR THE 
SCATTERING MATRIX 

In the past, quantum electronic transport in meso- 
scopic systems has been described successfully in terms 
of ensembles of single-electron S matrices (see, for in- 
stance, Refs. |,@,y,^|Jl6Hl|). In particular, it was 
shown that, in the absence of direct processes, quantum 
transport through classically chaotic cavities can be stud- 
ied in terms of the invariant measure for the S matrix, 
which is a precise mathematical formulation of the intu- 
itive notion of "equal-a-priori-probabilities" in S'-matrix 
space. For TRI systems, like the ones we are studying 
here, the invariant measure ip-also known as the Circular 
Orthogonal Ensemble (COE)K The COE for S matrices 
of dimensionality 2N, N being the number of open chan- 
nels supported by the two leads attached to the cavity, 
gives, for the ensemble averaged (indicated by .brackets 
(■ • ■)) spinless conductance and its varianceESlij 



(T) 



N 2 



var(T) 



2N + 1 



N(N + l) 2 



(2JV + iy (2N + 3)' 



(5.1) 



(5.2) 



respectively. The 1 in the denominator of Eq. (5.1) is 
the weak- localization correction (WLC). 

We construct here an ensemble of S matrices for the 
system consisting of a single electron and the PB, and, 
using the conductance formula obtained in the previous 
section, Eq. (4.34), we analyze the effect of the PB on 
the average and variance of the conductance. 

Within the model we have been discussing for the e- 
PB system, we postulat e m independent COE's for the 
S matrices S of Eq. (4.15). Ensemble averaging Eq. 



( 4.33 ) we obtain 



<T=l 

m 
<t=1 

N 2 



2N + 1 



2N+V 



(5-3) 



the same result as in Eq. (5.1), in the absence of the 
PB. Thus the model is not generic enough to decrease 
the WLC. For the variance we obtain 



•(T)] m = J2\p(°)} 2 var{tr[(ryr]} 



N(N + l) 2 



{2N+iy{2N + 3) ^ 



If one PB state exo has probability 1 and all the others 0, 
Y^=i [p( a )] 2 = I- In the other extreme case of equiprob- 
able PB states, Y^=i b( cr )] 2 = 1/m. The conductance 
thus fluctuates less than in the absence of the PB, as 
expectedB. 

Should the conjecture we made at the end of last sec- 
tion be true, and the result (4.34) be valid more gener- 
ally, i.e. for a generic e-PB S matrix, we could postulate 
a COE for the full 2miV-dimensional S matrix in the ND 
basis. We would then obtain, ensemble averaging Eq. 



(4.34) 



(T) m = 

fAVi/' ab 



if v 

z ab 



(5.5) 



From Ref. 12 we find 



5p V 



2mN+ 1 

and, since J2 U p u p B = 1, we finally obtain 

N 2 

(T), 



2N + — 



(5.6) 



Now we observe that, as m — > oo, the effect of the PB is 
to kill the WLC, as expected. 

We find the variance of T only in the simplified situa- 
tion described by Eq. (4.36) 



E E K 

_fi,i/—l ab 



(j.vi2 
b 



(5.7) 



The variance appearing on the right hand side of this last 
equation can be read form Eq. (5.2), with the replace- 
ment N => mN, to find 



[var(T)] f 



1 



1 ^,2 



\ m / V m J 



(5.8) 



For large m, conductance fluctuations are killed as well, 
as expected. 
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VI. CONCLUSIONS 



We have discussed the quantum electronic conduction 
problem in a mesoscopic system, allowing for the presence 
of a phase breaker (PB) the system can interact with. 
The PB can exist in m quantum-mechanical states. It 
could represent an impurity with internal degrees of free- 
dom, so that its state may change via the interaction with 
the electrons: for instance, it could be a magnetic impu- 
rity interacting with the electron spin. The PB might 
also represent the environment -for instance, the phonon 
field- whose state is allowed to change. 

We first studied the problem of a single electron in- 
teracting with the PB, in order to investigate the effect 
of the latter in the interference terms that are there in 
the absence of the PB. We described the static, as well 
as the dynamical scatterers (the PB), in terms of their 
scattering or S matrices: this makes the discussion very 
intuitive, quite general, and amenable to the application 
of random-matrix models that have been developed in 
the past. We found an S'-matrix formulation that is ca- 
pable of decribing the crossover from a purely coherent 
response to a classical, or incoherent one. 

We then set out to study the conduction problem -a 
many-electron problem. The e-e interaction is treated in 
an RPA approximation. We could make the very general 
statement (that so far, to our knowledge, was known only 
in the absence of a PB) that in the dc limit one can give 
an explicit expression for the current in terms of the po- 
tential difference applied between the two reservoirs, the 
full potential profile not being needed. That expression 
involves the conductivity tensor <r°^(r, r') which, in the 
absence of the PB, can be calculated in terms of single- 
particle Green's function and, eventually, single-particle 
transmission coefficients. In the presence of the PB, the 
e-e interaction induced by the PB has not allowed us 
to follow a similar path. That induced interaction is 
there even in the absence of the e-e Coulomb interac- 
tion, which, within RPA, was disposed of and replaced, 
in turn, by the solution of a self-consistent problem. 

We proposed a model for the e-PB interaction that can 
be diagonalized and transformed into m single-particle 
problems. The conduction problem splits into m single- 
particle problems as well, and the final result for the con- 



ductance, given in Eq. (4.34), is like the standard one 



without a PB, except that now an extra trace over the 
to PB states a ppear s. We should stress that, within the 
model, result (4.34) is not perturbative in the e-PB in- 
teraction strength. However, the model implies a single- 
electron-PB S matrix of a rather restricted form. For 
a "generic" single-electron-PB S matrix we do not even 
know under what approximations (that would imply dis- 
regarding the induced e-e interaction) the conduction 
problem can be reduce d to a single-electron one. The re- 
sult expressed in Eq. (4.34) is so intuitive, though, that 
we conjecture its validity under a more general e-PB S 
matrix, within some suitable approximation, as discussed 



at the end of section III around equation ( |3.17 ). 

It would be instructive to solve this same concLuction 
problem within the spirit of Landauer's approach!!! , jus t 
as in Ref. |[ and verify that one arrives at the result (4.34) 
under the special model used here and not in general. 
A suitable approximation for treating the more general 
problem might suggest itself in that approach. But this 
we have not succeeded to do so far. 

Finally, a random-matrix model was setup for the 
description of the e-PB system, with possible applica- 
tions to chaotic cavities: the effect of the PB on the 
conductance average and its fluctuations was analyzed. 
The limitations of the model became apparent in that 
study; in contrast, a generic e-PB S matrix was shown 
to give much more freedom in the description of the weak- 
localization correction and the conductance fluctuations. 
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